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Abstract 

We analyse the impact of explicit CP-violation in the Higgs sector of the Next-to-Minimal Su¬ 
persymmetric Standard Model (NMSSM) on its consistency with the Higgs boson data from 
the Large Hadron Collider (LHC). Through detailed scans of the parameter space of the com¬ 
plex NMSSM for certain fixed values of one of its CP-violating (CPV) phases, we obtain a 
large number of points corresponding to five phenomenologically relevant scenarios containing 
^ 125 GeV Higgs boson(s). We focus, in particular, on the scenarios where the visible peaks 
in the experimental samples can actually be explained by two nearly mass-degenerate neutral 
Higgs boson states. We find that some points corresponding to these scenarios give an overall 
slightly improved fit to the data, more so for non-zero values of the CPV phase, compared to 
the scenarios containing a single Higgs boson near 125 GeV. 
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1 Introduction 


The Higgs sector of the NMSSM [I] (see, e.g., las! for reviews) contains two additional neutral 
mass eigenstates besides the three of the Minimal Supersymmetric Standard Model (MSSM). This 
is due to the presence of a Higgs singlet superfield besides the two doublet superfields of the MSSM. 
When all the parameters in the Higgs and sfermion sectors of the NMSSM are real, one of these 
new Higgs states is a scalar and the other a pseudoscalar. Hence, in total three scalars, i?i, 2,35 
two pseudoscalars, Ai^ 2 ^ make up the neutral Higgs boson content of the model. This extended 
Higgs sector of the NMSSM boasts some unique phenomenological possibilities, which are either 
precluded or experimentally ruled out in the MSSM. For example, in the NMSSM either of the 
two lightest CP-even Higgs bosons. Hi or i72, can play the role of the ^ 125 GeV Standard Model 
(SM)-like Higgs boson, i^obs? observed at the LHC |ll[5l[6]. 

Of particular interest in the NMSSM is the possibility that the SM-like Higgs boson can obtain 
a large tree-level mass in a natural way, i.e., without requiring large radiative corrections from the 
supersymmetric sectors. This happens in a specific region of the parameter space, which we refer to 
as the natural NMSSM, where there is a significant singlet-doublet mixing and the Lfobs is typically 
H 2 . This scenario was used to explain [7] the enhancement in the Lfobs ^ 77 channel in the early 
LHC data. However, when the singlet-doublet mixing is too large, the properties of H 2 can deviate 
appreciably from an exact SM-like behaviour, resulting in a reduction of its fermionic partial decay 
widths. An alternative possibility in a very similar parameter space region is that of both Hi 
and H 2 simultaneously having masses near 125 GeV UM- In that case, the observed excess at 
the LHC could actually be due to a superposition of these two states, when their individual signal 
peaks cannot be resolved separately. One of these two Higgs bosons, typically iLi, is the singlet-like 
neutral state. Moreover, in m it was noted that the lighter of the two pseudoscalars, Ai, when it 
is singlet-like, could also be nearly mass-degenerate with a SM-like Hi near 125 GeV, instead of or 
even along with the H 2 . However, such a pseudoscalar can only contribute visibly to the measured 
signal strength near 125 GeV if it is produced in association with a bb pair. 

One of the most important yet unresolved issues in particle physics is that of the observed 
matter-antimatter asymmetry in the universe. A plausible explanation for this asymmetry is elec- 
troweak (EW) baryogenesis [TT]. The necessary conditions for successful EW baryogenesis include 
the following [12]: (1) baryon number violation, (2) CP-violation and (3) departure from equilib¬ 
rium at the critical temperature of the EW symmetry breaking (EWSB) phase transition, implying 
that it is strongly first order. In the SM, a strongly first order EW phase transition is not possible 
given the measured mass of the Higgs boson at the LHC. Besides, the only source of CP-violation 
in the SM, the Cabibbo-Kobayashi-Maskawa matrix, is insufficient. Therefore, beyond the SM, a 
variety of sources of CP-violation have been proposed in the literature (for a review, see [I3| and 
references therein). In the context of supersymmetry (SUSY), a strongly first order phase transition 
is possible in the MSSM only if the lightest stop has a mass below that of the top quark. This 
possibility has now been ruled out by SUSY searches at the LHC [T4|. Also, the MSSM Higgs 
sector does not violate CP at the tree level but does so only at higher orders [I5l[l6l[17|. The 
CPV phases, transmitted radiatively to the Higgs sector via couplings to the sfermions, are tightly 
constrained by the measurements of fermion electric dipole moments (EDMs) [TH|. However, these 
EDM constraints can be relaxed under certain conditions [MiiiniEoi. 

The NMSSM has been shown to accommodate a strongly first order EW phase transition with¬ 
out a light stop [21]. Additionally, in this model, CP-violation can be invoked explicitly in the Higgs 
sector even at the tree level by assuming the Higgs self-couplings, A and to be complex. Beyond 
the Born approximation, the phase of the SUSY-breaking Higgs-sfermion-sfermion couplings, Ap 
where / denotes a SM fermion, is also induced in the Higgs sector, as in the MSSM. In the presence 
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of the associated complex phases, the five neutral Higgs bosons are CP-indefinite states, due to 
the mixing between the scalar and pseudoscalar interaction eigenstates. CPV phases can therefore 
influence the phenomenology of the NMSSM Higgs bosons by, e.g., modifying their mass spectrum 
as well as their production and decay rates [22], similarly to the MSSM [23|. The impact of these 
phases in the complex NMSSM (cNMSSM), i.e., the CPV NMSSM, on the necessary conditions 
for successful EW phase transition was also studied some time ago m- The consistency of sce¬ 
narios yielding the correct baryon asymmetry with the LHC Higgs boson data still remains to be 
studied in depth, though. However, even leaving aside these considerations, the possibly distinct 
phenomenological scenarios that the cNMSSM can yield make it a particularly interesting model 
for exploration at the Run-H of the LHC. 

The cNMSSM has therefore been the subject of several studies recently and, in particular, 
some important theoretical developments have been made in the model. The dominant 1-loop 
corrections to the neutral Higgs sector from the (s)quark and gauge sectors were studied in [251126] . 
in the renormalisation group equations-improved effective potential approach. The corrections from 
the gaugino sector were included in m and, more inclusively, recently in [2H|. In the Feynman 
diagrammatic approach, the complete 1-loop Higgs mass matrix was derived in [29| and the 0{atas) 
contributions to it were calculated in [301 . As far as the phenomenology of the Higgs bosons in the 
cNMSSM is concerned, the consistency of several CPV scenarios with the early results on the Lfobs 
from the LHC data was studied in detail in pniES]. Another distinct phenomenological scenario, 
possible only for non-zero CPV phases, has also been studied in m- 

The CMS and ATLAS collaborations have recently updated their measurements of the Lfobs 
signal rates in the and bb channels [SUES]. The fact that these rates also tend to favour 

a SM-like Lfobs is increasingly jeopardising the above mentioned natural NMSSM scenario with 
large singlet-doublet mixing but only with one Higgs boson, either Hi or iL 2 , around 125 GeV. This 
makes the scenario with both Hi and H 2 contributing to the observed ^ 125 GeV signal all the more 
important, since it may potentially satisfy better the current Higgs boson data while still leaving 
plenty of room for new physics. In case of the cNMSSM, since the five neutral Higgs bosons are 
CP-mixed states, the scenario with mass-degenerate Hi and H 2 can entail both the corresponding 
possibilities in the real NMSSM (rNMSSM), i.e., mass-degenerate Hi, H 2 or Hi/H 2 ^ Ai. 

In this study we therefore analyse and compare the prospects for scenarios with two mass- 
degenerate Higgs bosons against those with a single Higgs boson near 125 GeV in the Za-invariant 
cNMSSM. We perform scans of the relevant parameter space [8] of the model using the public 
program NMSSMCALC [33| to search for all possible ^ 125 GeV Higgs boson scenarios, with the 
CPV phase of the coupling k set to five different values, including 0° - the rNMSSM limit - in each 
case. The condition for mass-degeneracy between two Higgs bosons is imposed by requiring them to 
lie within 2.5 GeV of each other, which is consistent with the current mass resolution of the LHC [M] . 
taking into account the uncertainties in the theoretical mass prediction. We then use fits to the 
Higgs boson data from the LHC Run-I, both with = 7TeV and = 8TeV, as well as from 
the Tevatron, performed using the program HiggsSignals [35], as the sole criterion for comparing 
the present likelihood of each of these scenarios. We also discuss how these mass-degenerate Higgs 
bosons can be identified at the LHC based on the signal rate double ratios introduced in [36] . 

The article is organised as follows. In the next section we will briefly revisit the Higgs sector 
of the cNMSSM. In section [3] we will provide details of our numerical scans and our procedure for 
fitting the model predictions for the Higgs boson(s) to the LHC data. In section 0] we will discuss 
the results of our analysis and in section [5] we will present our conclusions. 
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2 The Higgs sector of the cNMSSM 

The NMSSM contains a singlet Higgs superfield, S, besides the two SU{2)l doublet superfields, 


Hu = 


H+ 

H?, 



( 1 ) 


of the MSSM. The superpotential of the NMSSM is written as 

h^NMSSM = MSSM Yukawa terms + XSHuHd + — , (2) 

o 

where A and k are dimensionless Yukawa couplings. This superpotential is scale invariant, since the 
term fjiHuHd appearing in the MSSM superpotential has been removed by imposing a discrete Z 3 
symmetry. In this model, an effective /i-term, /ieff = As, is instead generated when the singlet field 
acquires a vacuum expectation value (VEV), s, which is naturally of the order of the SUSY-breaking 
scale. 

The tree-level Higgs potential of the NMSSM, obtained from the superpotential in eq. ([2]), is 
written in terms of the neutral scalar components of the Higgs superfields, and 5, as 

Ho = 

+ + |A5|2) 

+X {K + - \< - +1 

+m||5|2 + {XAx {H+HJ - hM) S + + h.c.) , (3) 

where , with gi and ^2 being the U{1)y and SU{2)l gauge couplings, respectively, and 

Ax and are the soft SUSY-breaking Higgs trilinear couplings. The scalar fields and S 

are developed around their respective VEVs, Vu^ Vd and 5 , as 


s 



{vd + HdR + iHdi) 


h: 


A6 


Hi 


^^2 HuR T 


ul) 


gi<p 

71 


(s + Sr + iSj). 


( 4 ) 


The Higgs coupling parameters appearing in the potential in eq. (jSl) can very well be complex, 
implying A = |A|e*'^^, k = |K|e*'^'', Ax = \A\\e^'^^^ and . As a result, the Vq) 

evaluated at the vacuum, contains the phase combinations 


4>x = 4>\ + ^ I 4 >k = 4>k + ^ 4>x + 4 >Ax (f)'^ + (f)A^ ■ 


( 5 ) 


For correct EWSB, the Higgs potential should have a minimum at non-vanishing Vd and s, 
which is ensured by requiring 

= 0 for $ = HdR, HuR, Sr, Hdi, Hui, Si. 



4 


( 6 ) 



Through the above minimisation conditions the phase combinations (j)'^ + (j)Ax 0'^ + (j)A^ can be 
determined up to a twofold ambiguity by (j)'^ — Thus, (j)'^ — (j)'^ is the only physical CP phase 
appearing in the NMSSM Higgs sector at the tree level. Also, using these conditions, the soft mass 
parameters traded for the corresponding Higgs field VEVs. 

The neutral Higgs mass matrix is obtained by taking the second derivative of the Vq evaluated 
at the vacuum. This 5x5 matrix, A^q, in the HuR^ Sr^ Hr Sj) basis, from which the 

massless Nambu-Goldstone mode has been rotated away, can be diagonalised using an orthogonal 
matrix, O, as O^M.qO = diag(m|^^ ^1^5)• This yields the physical tree-level masses 

corresponding to the five mass eigenstates, 

{Hi, H2, Hs, Hi, HAi = Oai {HdR, Hur, Sr, Hi, Si)J , (7) 

such that The elements, of fhc mixing matrix then govern 

the couplings of the Higgs bosons to all the particles in the model. 

The tree-level Higgs mass matrix is subject to higher order corrections from the SM fermions, 
from the gauge and chargino/neutralino sectors and the Higgs sector itself, as well as from the 
sfermion sector, in the case of which they are dominated by the stop contributions. Upon the 
inclusion of these corrections, AAd^, the Higgs mass matrix gets modified, so that 

= + AM^ . ( 8 ) 

Explicit expressions for A4q as well as for AA4^ can be found in [29l EZl EE]- Thus, beyond the 
Born approximation, the CPV phases of the gaugino mass parameters, Mi^ 2 , and of A^ are also 
radiatively induced in the Higgs sector of the NMSSM. 

Therefore, when studying the phenomenology of the Higgs bosons, one needs to take into 
account also the parameters from the other sectors of the model. However, the most general 
NMSSM contains more than 130 parameters at the EW scale. Assuming the matrices for the 
sfermion masses and for the trilinear scalar couplings to be diagonal considerably reduces the 
number of free parameters. One can further exploit the fact, mentioned above, that the corrections 
to the Higgs boson masses from the sfermions are largely dominated by the stop sector. For our 
numerical analysis in the following sections, we will thus impose the following supergravity-inspired 
universality conditions on the model parameters at the EW scale: 

Mo = Mg^ 2,3 = ^Ui,2,3 = ^Di,2,3 = ^Li^2,3 = ^Ei,2,3 ^ 

^ 1/2 = 2Mi = M 2 — gTfs , (9) 

Aq = A^ = = Af , 

where and are the squared soft masses of the sfermions, 

Tfi,2,3 those of the gauginos and A^^ ~ the soft trilinear couplings. Altogether, the input parameters 
of the cNMSSM then include 

Mq , IM 1 / 2 I, |Ao|, tan^ (= Vu/vd ), |A|, \r\ , /ieff, |^a| , , O 1/2 , 6^/, 0'^ , 

where 61/2 and 9^ are the phases of the unified parameters Mx /2 and Aq, respectively. 

3 Numerical analysis 

As noted in the Introduction, non-zero CPV phases can modify appreciably the masses and decay 
widths of the neutral Higgs bosons compared to the CP-conserving case for a given set of the 
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remaining free parameters. In the case of the i^obs candidate in the model, whether Hi or H 2 
or even the CPV phases are thus strongly constrained by the LHC mass and signal rate 
measurements. This was analysed in detail in [ 22 ] , where the scenarios with mass-degenerate Higgs 
bosons were, however, not taken into account. In the present study we thus test whether the said 
modifications in the Higgs boson properties with non-zero values of the phase 0'^ (by which we 
imply 0^, which is the actual physically meaningful phase, since ip can be absorbed into (/>(. by a 
field re-definition) can lead to a relatively improved consistency with the experimental data. 

The reason for choosing (j)'^ as the only variable phase, while setting ^ 1 / 2 , 0? and to 0°, is 
that it is virtually unconstrained by the measurements of fermionic EDMs [26l[29]. Furthermore, 
our aim here is to analyse the scenarios with a generic CPV phase and compare them with the 
rNMSSM limit rather than measuring the effect of any of the individual phases. Note however that, 
since only the difference (j)'^ — (/)(. enters the Higgs mass matrix at the tree level, the impact of a 
variation in 0'^ is also quantified by that due to the variation in 0'^ at this level. At higher orders 
though, a variation in has an impact on the sfermion and neutralino/chargino sectors which is 
independent of 

In our numerical analysis, we used the program NMSSMCALC-vl.03 [33| for computing the 
Higgs boson mass spectrum and decay branching ratios (BRs) for a given model input point. The 
public distribution of NMSSMCALC contains two separate packages, one for the rNMSSM only and 
the other for the cNMSSM. Some supersymmetric corrections to the Higgs boson decay widths are 
currently only available in the rNMSSM and hence are not included in the cNMSSM package. For 
consistency among our rNMSSM and cNMSSM results, we therefore set = 0° in the cNMSSM 
package for the rNMSSM case instead of using the dedicated rNMSSM package. Furthermore, 
using the cNMSSM code also for the rNMSSM limit makes it convenient to draw a one-on-one 
correspondence between the = 0° case and each of the > 0° cases in a given scenario. This 
is because in the cNMSSM package, even in the rNMSSM limit, the five neutral Higgs bosons are 
ordered by their masses and not separated on the basis of their CP-identities. Thus, the scenario 
with mass-degenerate iLi, iL 2 , which we will henceforth refer to as the iLobs — Hi + H 2 scenario, 
takes into account both the ^ 125 GeV iLi, H 2 as well as the ^ 125 GeV iLi, Ai solutions of 
the rNMSSM without distinguishing between them. If one, conversely, uses the rNMSSM package, 
these two scenarios ought to be considered separately. The same is true also for the iLobs = ^2 + ^3 
scenario, wherein iL 2 , H 3 are mass-degenerate. 

The program NMSSMCALC allows one the option to include only the complete 1 -loop contribu¬ 
tions in the Higgs mass matrix or to add also the 2-loop O(atas) corrections to it. For our analysis, 
in order for better theoretical precision, we evaluated the Higgs boson masses at the 2-loop level. In 
the NMSSMCALC input, one also needs to choose between the modified dimensional regularisation 
(DR) and on-shell renormalisation schemes for calculating contributions from the top/stop sector 
in the program. We opted for the DR scheme for each scenario. Note though that further inclusion 
of 0{ai)as)^ 0{at + 0^5 + o^r)^ and the recently calculated NMSSM-specific 0{ax -h 2-loop 
corrections m in NMSSMCALC may have a non-negligible impact on the Higgs boson masses and 
observables [3H|. We, however, maintain that such contributions will only result in a slight shifting 
of the parameter configurations yielding solutions of our interest here, but our overall results and 
conclusions should still remain valid. 

We performed six sets of scans of the cNMSSM parameter space by linking NMSSMCALC 
with the MultiNest-v2.18 [39] package. MultiNest performs a multimodal sampling of a theoretical 
model’s parameter space based on Bayesian evidence estimation. However, we use this package 
not for drawing Bayesian inferences about the various NMSSM scenarios considered but simply to 
avoid a completely random sampling of the 9-dimensional model parameter space. In the program, 
we therefore defined a Gaussian likelihood function for the iLobs ^ given scan, assuming the 
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Parameter 

Natural NMSSM range 


Parameter 

Low-A-NMSSM range 

Mo (GeV) 

200 - 2000 


Mo (GeV) 

200 - 4000 

Mi /2 (GeV) 

100 -1000 


Mi /2 (GeV) 

100 - 2000 

^0 (GeV) 

-3000 - 0 


^0 (GeV) 

-7000 - 0 

tan/3 

1-8 


tan (3 

5-45 

A 

0.4 - 0.7 


A 

0.001 - 0.4 

K 

0.3 - 0.6 


K 

0.001 - 0.3 

Meff (GeV) 

100 - 300 


l^es (GeV) 

100 - 2000 

AIa (GeV) 

-1000 - 1000 


(GeV) 

-1000 - 4000 

AIk (GeV) 

-1000 - 1000 


(GeV) 

-4000 - 1000 


(a) (b) 


Table 1: Ranges of the NMSSM parameters scanned, with fixed 0^, for (a) each i^obs scenario in 
the natural NMSSM and (b) the low-A-NMSSM scenario. 


experimental measurement of its mass to be 125 GeV and allowing upto ±2 GeV error in its model 
prediction. We set the enlargement factor reduction parameter to 0.3 and the evidence tolerance 
factor to a rather small value of 0.2, so that while the package sampled more concentratedly near the 
central mass value, a sufficiently large number of points were collected before the scan converged. 
In each of the first two sets of scans we required Hi to be the i^obs- Iii fhe third set we imposed 
this requirement of consistency with the i^obs mass on i72, in the fourth on in the fifth on both 
i^i, H 2 and in the sixth on both i72, H^. Each of the six sets further contained five separate scans 
corresponding to = 0° , 3° , 10° , 30° and 60°. 

The scanned ranges of the nine free parameters (after fixing the phases) of the natural NMSSM, 
which are uniform across all its five scenarios considered, are given in tab. [I](a). Only large values 
of A and k. are used in this model (with the upper cut-off on them imposed to avoid the Landau 
pole). Since large radiative corrections from SUSY sectors are not necessary in the natural limit 
of the NMSSM, the parameters Mq, M 1/2 and Aq are not required to take too large values. Note 
that while Aq can in principle be both positive and negative, with a slightly different impact on 
the physical mass of the SM-like Higgs boson for an identical set of other input parameters in each 
case, we restricted the scans to its negative values only, in order to increase the scanning efficiency. 

In the remaining sixth scan, we considered the complementary parameter space of the NMSSM, 
with A and k, kept to relatively smaller (and tan (3 to larger) values, so as to prevent too large a 
singlet-doublet mixing. In fact, for A, ^ 0, when the singlet sector gets effectively decoupled, 
iLi, which is by default identified with the iLobs? has properties very identical to the lightest Higgs 
boson of the MSSM. Since Hi in such a case does not obtain a maximal tree-level mass that is 
possible in the most general model, large radiative corrections are needed from the SUSY sector. 
Hence we used slightly extended ranges of the remaining parameters, which are given in tab. [Il^b). 
This scenario, which we refer to as the low-A-NMSSM scenario henceforth, has been included in 
our analysis in order to compare the inferences made for the natural NMSSM with an approximate 
MSSM limit of the model. 

Once the scans had completed, we filtered the points obtained with each by further imposing 
123 GeV < GeV. Note that in the iLobs — Hi -\- H 2 and iLobs = ^2 + ^3 scenarios, 

this condition was imposed on iL 2 , since in both these scenarios it is typically the Higgs boson 
with SM-like couplings. The total number of points, Vtotab remaining after this filter is given in 


7 








Scenario 

low-A 

^obs — Hi 

Hohs — H2 

^obs — H^ 

^obs — Hi + H2 

Hohs — H2 + 





w = o° 



Atotal 

17786 

15675 

15072 

14431 

26045 

23736 

Nub 

17722 

13691 

2904 

965 

11878 

2819 




(p 

•k = 3° 




17829 

15775 

15026 

14806 

27199 

25684 

-/Vhb 

17782 

13885 

3235 

2391 

11863 

1659 




</>/ 

. = 10° 




17847 

15784 

15080 

14810 

26735 

28348 


17786 

13866 

2411 

2495 

12607 

3369 




(p, 

0 

0 

CO 

II 




17810 

16256 

15037 

14671 

31719 

28685 

-^HB 

17743 

14725 

247 

276 

13503 

2012 




(pi 

. = 60° 




17810 

0 

14996 

14438 

0 

30412 

-/Vhb 

17743 

0 

247 

2 

0 

242 


Table 2: Number of scanned points remaining after imposing the mass constraint on i^obs those 
passing the HiggsBounds test, for each scenario studied. See text for details. 


tab. [2] for each scenario considered. All these points were then tested for consistency with the 
LEP and LHC exclusion limits on the other, non-SM-like, Higgs bosons of the model, using the 
package HiggsBounds v4.2.0 SOI. The points passing the HiggsBounds test were retained as the 
‘good points’ for further analysis, and their number, denoted by A^hb, for each scenario is also given 
in tab. [2l We point out for later reference that in each of the two Lfobs = Hi scenarios as well as 
in the Lfobs — Hi-\- H 2 scenario, the number of surviving good points (where they are available) is 
very identical across all input values of 0^, implying mutually fairly consistent sample sizes. 

Next we carried out fits to the Lfobs data for the good points using the public code HiggsSignals 
vl.3.2 [35]. For obtaining these fits, HiggsSignals requires, along with the masses and BRs of each 
Higgs boson, the square of its normalised effective couplings, {gUiX/dhsmxY^ ^ given SM 
particle pair X, with hsu being the SM Higgs boson with the same mass as the Hi. Note that 
when X is a pair of fermions, there is a scalar as well as a pseudoscalar normalised coupling for 
each Hi^ both of which need to be passed separately to HiggsSignals. All these are then used to 
calculate the normalised cross sections. 


X ^ ^{pp Hj^ X) 

“ a{pp -> hsM -> X) ’ 


( 10 ) 


corresponding to a given decay channel, X, in an approximate way. The NMSSMCALC version we 
used did not provide the normalised Higgs boson couplings as an output. We therefore modified 
the code to obtain these couplings for adding them as a dedicated block in the SLHA input file for 
HiggsSignals. 

The program HiggsSignals compares the computed for each Hi with the experimentally 
measured ones, /i^p, for wide ranges of input Higgs boson masses in a variety of its production and 
decay channels at the LHC and the Tevatron. We used only the ‘peak-centred’ method and the 
‘latestresults’ observable set in the program, with the assignment range variable A set to the default 
value of 1. It thus performed a fit to a total of 81 Higgs boson peak observables (77 from signal 
strength and 4 from mass measurements), from the CMS, ATLAS, CDF and D0 collaborations. 



















































for a given model point. We assumed a Gaussian theoretical uncertainty of 2GeV in the masses 
of the three lightest neutral Higgs bosons of the model. The default values of the uncertainties in 
the Higgs boson production cross sections as well as BRs were retained. Further details about the 
fitting procedure can be found in the manual [35] of the package. The main output of HiggsSignals 
contains the total and the p-value from the fit, given the number of statistical degrees of freedom, 
for each model point. Since the aim of this study is a comparison of various i^obs scenarios rather 
than the overall goodness of fit for each, we will quantify our results only in terms of the and 
ignore the p-value. 

As an observable indication of the presence of more than one Higgs bosons near 125 GeV, the 
double ratios 


Di^ 


^vbf(77)/-^^^(77) 

KsAbb)/RU^b) 


D2 = 


R'YBFi'y'y) / ^ggi'y^) 

KMWW)/R^g{WW) 


D3 = 


R'^BFiWW)/R^g{WW) 

Rh^^{hb)/Rhg{bb) 


,( 11 ) 


were proposed in [30] • Each of these ratios should be unity if the -ffobs constitutes of only a single 
Higgs boson, while the contribution of two (or more) Higgs bosons to the -Hobs signal could result 
in a deviation of these ratios from 1. In the above expressions, i?y(X) = Ry^X) + where 

Hi and Hj are the two mass-degenerate Higgs bosons in a given scenario and the subscripts VBF 
and gg imply the vector boson fusion and the gluon fusion production modes, respectively. 
for each Hi is defined as 


^ T{Hi ^ y) BR(H, ^ X) 
^ ^ ^ r(hsM ^ E) BR(hsM ^ E) 




( 12 ) 


with Y being the given production mode and, in the last equality, = ^{Hi Y (E))/r(hsM 

X{Y)), the normalised partial decay width of Hi into the X (Y) pair])] and E/jg^ are the total 
decay widths of Hi and /isM, respectively. 

We also evaluated the ratios T^i, D 2 and for the points which give reasonably good fits to 
the data (to be defined later) in the scenarios with two mass-degenerate Higgs bosons. For this 
purpose, for each Hi was calculated by fixing F/^g^^ in eq. (|T^ to 4.105 x 10“^ GeV, which 

is the value given by the program HDECAY [H] for a 125 GeV /isM- A change of ±2 GeV in the 
mass of /isM has only a marginal affect on this width, which we ignore. For calculating the 
with HDECAY, care was taken that all the partial decay widths of hsu were evaluated at the same 
perturbative order as that implemented in NMSSMCALC for computing F/f.. Moreover, Cy^ is 
simply the squared normalised coupling of Hi to a vector boson, V, pair for the VBF production 
mode and to a gluon pair for the gg mode. Similarly, implies the HiVV and Hi^^ normalised 
couplings squared, respectively, for X = WW and 77 . All these couplings are thus the same ones 
obtained from NMSSMCALC for passing to HiggsSignals. In the case of X = 56, though, there 
is a scalar and pseudoscalar coupling for each Hi^ as noted above. For this reason, G^^"’s were 
calculated using the actual T{Hi bb) from the NMSSMCALC output for a given model point 
and the F(/ism ^ bb) obtained from HDECAY for — 125 GeV. 


4 Results and discussion 

In fig. [I] we show the total obtained for the points from our scans for the various iLobs scenarios 
considered. The green points in the figure correspond to = 0°, violet to = 3°, blue to 

^Note that eq. m assumes that the /isM-normalised production cross sections for the Y = VBF and gg processes 
can be approximated by the normalised partial decay widths of Hi in the VV and gg decay channels, respectively. 
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= 10°, red to (p,^ = 30° and cyan to = 60°. For the scenarios in which only one of the three 
lightest neutral Higgs bosons is assumed to be the i^obs? we have made sure that the difference 
between the mass of i^obs that of each additional Higgs boson nearest to it is always larger 
than 2.5 GeV. The lower cut-off in each panel, in this figure and in those that follow, varies 

depending on the minimum value obtained in the corresponding scenario. The upper cut-off in 
for each scenario is chosen so as to include as many points in the corresponding figures as possible 
without the getting more than 10 units larger than the minimum obtained in that scenario 
(given that there are 9 statistical degrees of freedom). 

Fig. lHa) corresponds to the low-A-NMSSM scenario. One notices in the figure that for = 0° 
the Xmin des very close to 70, and is thus almost identical to Xmin ~ 69.96 that is given by 
HiggsSignals for a SM Higgs boson at a mass of 125.1 GeV, with the same settings as used by us. 
The input parameters (with the exception of Mq, Mi /2 Aq^ which can be adjusted with much 
more freedom) and the masses of the three lightest Higgs bosons are given in tab. [3l The negligibly 
small difference in the Xmin value obtained for the /ism and for the CP-conserving low-A-NMSSM 
results from the fact that A for the corresponding point in the latter is non-vanishing, as seen in 
the table, so that the singlet sector is not completely decoupled and an exact MSSM-limit is not 
reached. One can notice in the figure and the table a slightly lower value of Xmin obtained for the 
sets of points corresponding to non-zero values. However, A for all these points is even larger 
than in the CP-conserving limit. Note also that, for all 0^, most of the points give < 1. 

In fig. [D^b), which corresponds to the i^obs = Hi scenario in the natural NMSSM, we see that 
there is a large concentration of points above a x^ value which is very similar to the Xmin ^be 

adjacent fig. [U^a), for each corresponding For non-zero though, one also sees a few scattered 
points with x^ lower than that for any of the points in the high concentration region. The overall 
lowest x^ lies very close to 68, for = 30°, with the mass of i^obs foi* fhe corresponding point lying 
at 124.5 GeV. However, according to tab. [3l the mass of H 2 for this point is within 3 GeV of that 
of Hi. It is therefore very likely that the relatively better fit for this particular point is a result of 
the assignment of H 2 instead of or along with Hi to some of the observables, especially when their 
experimental mass resolution is relatively poor. This possibility, which implies that our assumption 
of two Higgs bosons being individually irresolvable if their masses lie within 2.5 GeV of each other 
is rather robust, will be discussed further later. For = 60° none of the points obtained in the 
scan for this scenario had Hi heavier than 123 GeV. 

In the i^obs = H 2 scenario, a much smaller number of points was passed by HiggsBounds 
compared to the i^obs = Hi scenario, as seen in Fig. [D^c), but the equally low for most 

here, including 0°. Only for — 60°, while plenty of points with mH 2 ~ 125 GeV were obtained 
in the scan, the x^ foi* them is never low enough to appear in the figure. Once again, in tab. [3] one 
can see that, for the points giving the lowest x^ foi* each in this scenario. Hi always lies within 
3-4 GeV of H 2 . Hence the slightly better fit for this point is again made possible by a contribution 
of Hi to some search channels. In fig. [D^d) for the i^obs = H^ scenario, although very few points 
with Ax^ < 10 appear in this scenario compared to the ones above, the Xmin similar, except 

for the 0/^ = 0° case, when it has a fairly high value of around 77. 

In fig. [U^e) is shown the total x^ for the i^obs = Hi + H 2 scenario against the H 2 mass. One can 
observe quite a few similarities between this figure and the fig. [Ifb) seen above (for the i^obs = Hi 
scenario). There is once again a large concentration of points with x^ ^ 69 for all except 60°, 
and also many scattered points below it. Importantly though, there are many points in this scenario 
which give a x^ lower than 68, which is the overall lowest value observed for any other scenario 
here. Most of these points, including the one with the overall lowest x^ of 65, correspond to 
0/^ = 10°, although some points for other can also be noticed. In fig. [T]^f) one sees a Xmin ^f 
68 for the i^obs — H 2 A H^ scenario also but very few points with x^ < 71? in contrast with the 
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Figure 1: Total a function of: (a, b) when only ffi is assumed to be the i^obs? foi* fwo 

different sets of scanned ranges of the parameter space (see text for details); (c) when only H 2 
is the i^obs; (d) m/fg when only is the i^obs; {^) '^H 2 when both Hi and H 2 lie near 125 GeV 

and (f) 171^2 when both H 2 and lie near 125 GeV. 
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Scenario 

low-A 

^obs — Hi 

^obs — H 2 

^obs — 

^obs — Hi + H 2 

Hobs — H 2 + H 3 




A 

c = 0° 



Xmin 

70.1 

69.5 

68.5 

76.9 

65.9 

69.8 

A 

0.046 

0.582 

0.653 

0.48 

0.597 

0.597 

Hi 

0.213 

0.43 

0.511 

0.305 

0.302 

0.327 

tan (3 

17.65 

1.66 

3.6 

6.98 

2.39 

2.07 


853.6 

226.8 

609.7 

680.7 

540.0 

179.3 

A. 

-2352 

-741.4 

—666.0 

14.05 

-479.3 

-3.95 

Meff 

130.0 

281.5 

243.7 

102.6 

285.2 

112.3 

rriHi 

125.4 

125.3 

122.1 

66.8 

123.3 

115.1 

mH2 

162.8 

142.1 

125.1 

121.0 

125.5 

125.1 

rriHs 

1828 

510.6 

618.5 

125.7 

730.0 

126.6 





c = 3° 



Xmin 

69.7 

69.2 

68.1 

68.2 

66.0 

68.1 

A 

0.184 

0.639 

0.588 

0.662 

0.631 

0.636 

Hi 

0.291 

0.523 

0.39 

0.349 

0.373 

0.318 

tan (3 

29.6 

1.81 

2.61 

4.24 

1.61 

6.45 

Ax 

2175 

162.5 

459.6 

425.6 

222.0 

848.6 

A. 

-236.7 

-595.1 

-597.6 

-12.03 

345.2 

-19.4 

Meff 

177.9 

218.8 

260.5 

110.1 

196.4 

127.4 

rriHi 

125.1 

125.3 

122.5 

97.2 

123.4 

105.1 

mH2 

444.9 

141.7 

125.8 

122.3 

125.2 

125.0 

rriHs 

496.1 

405.5 

563.6 

126.0 

366.3 

127.2 




(pK 

= 10° 



Xmin 

69.7 

68.8 

69.0 

69.4 

65.1 

68.1 

A 

0.138 

0.68 

056 

0.692 

0.688 

0.585 

K 

0.219 

0.409 

0.345 

0.338 

0.361 

0.306 

tan/3 

16.7 

1.85 

1.91 

4.88 

1.98 

7.55 


1379 

291.6 

347.8 

557.0 

390.8 

972.6 

A. 

-623.8 

-476.1 

-567.8 

12.7 

-435.1 

-30.62 

Meff 

133.5 

251.0 

266.9 

124.3 

254.0 

136.7 

rriHi 

125.0 

125.3 

120.3 

106.4 

123.6 

118.7 

mH2 

212.2 

140.5 

124.5 

111.6 

126.0 

126.1 

rriHs 

631.6 

482.5 

541.8 

125.6 

440.1 

127.4 




(pK 

= 30° 



Xmin 

69.7 

68.1 

68.6 

70.4 

65.6 

70.2 

A 

0.136 

0.648 

0.679 

0.537 

0.624 

0.481 

K 

0.219 

0.319 

0.586 

0.303 

0.388 

0.311 

tan/3 

29.4 

2.2 

2.13 

6.55 

2.10 

7.67 

^A 

3515 

570.1 

295.0 

702.2 

345.7 

796.5 

A. 

-781.0 

-398.4 

-590.7 

7.07 

330.5 

-23.22 

Meff 

170.8 

288.5 

227.8 

112.6 

209.1 

110.0 

niHi 

125.1 

124.5 

123.1 

86.5 

121.6 

107.1 

mH2 

234.3 

127.4 

126.1 

116.8 

123.8 

124.7 

mHz 

857.7 

462.4 

507.8 

124.3 

405.8 

125.8 


Table 3: Input parameters and Higgs boson masses corresponding to the points giving the lowest 
for all cases in each of the i^obs scenarios considered. 
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i^obs = Hi and i^obs — Hi-\- H 2 scenarios but similarly to the i^obs = ^2 and i^obs = ^3 scenarios. 

From the above discussion, it is clear that certain points, or parameter space configurations, in 
the i^obs = H 1 + H 2 scenario give the best fit to the current experimental Higgs boson data. A global 
Xmin’ lowest value across all scenarios examined here, of around 65 has been observed 

for (j),^ — 10° in this scenario, with some points corresponding to other values of also lying within 
1 unit of this None of the points obtained for the other scenarios gives a lying even within 
3 units of this global minimum, despite the number of sampled points for the i^obs = Hi scenario 
being typically larger. The reason for a better fit for some points with two nearly degenerate Higgs 
bosons becomes apparent by looking at the detailed output of HiggsSignals. In the peak-centred 
method, HiggsSignals assigns to a given observable the Higgs boson with a mass closest to the 
measured mass provided by the experiment. This mass measurement currently ranges between 
124.7 GeV to 126.0 GeV. Thus, when a single Higgs boson is assigned to all the observables, the 
contribution is large from the observables for which the measured mass lies away from the mass 
of the assigned Higgs boson, and the experimental mass resolution is good. On the other hand, 
when two Higgs bosons lie close to each other, the one assigned to a given observable is the one 
for which the difference of the predicted mass from the experimental value is the smallest, so that 
the contribution from this observable is minimal. This is as long as the mass of the other Higgs 
boson nearby lies outside the experimental mass resolution, otherwise HiggsSignals automatically 
assigns both the Higgs bosons to an individual observable if it improves the fit. 

Some caveats are in order here though. A Ay^ 2 ^ 3 is statistically quite insignificant for drawing 
any concrete inferences about the considered scenarios, since the total number of observables and 
statistical degrees of freedom is quite large. At the same time, the number of points giving Ay^ < 3 
is also fairly small. Moreover, no other experimental constraints have been imposed in our analysis, 
since the publicly available tools for testing these are so far not compatible with the cNMSSM. 
It is thus possible that many of the interesting points may have already been ruled out by such 
constraints. However, the aim of this study is not to disregard one scenario in favour of another, but 
to simply show that, given the current experimental data, the scenario with two mass-degenerate 
Higgs bosons in the NMSSM provides as good, if not better, a fit as the scenarios with a single 
Higgs boson near 125 GeV. This alternative possibility even points towards a source of CP-violation 
beyond the SM and, therefore, warrants more dedicated analyses as well experimental probes. In 
the following we discuss some other interesting aspects of this scenario. 

In the left, middle and right panels of fig. [2] we show the ratios T^i, D 2 and 7^3, respectively, as 
functions of the mass difference, mH 2 ~ foi" various values in the i^obs — Hi-\- H2 scenario. 
The heat map corresponds to the total y^ obtained for the points shown in each panel. This y^ has 
a uniform upper cut-off of 71 across all panels, as in fig. [D^e), but its lower cut-off varies according 
to the minimum obtained for the (j),^ case that a given panel corresponds to. According to fig. Ha), 
for = 0° the three ratios remain largely close to unity, but deviations up to 15 — 20% can be seen 
for some points. D 2 ^ the ratio dependent on only the bosonic signal strengths, only gets smaller 
than 1 for some points and its maximum observed deviation is lower than that of Di and 7^3, each 
of which can be both above or below unity. Importantly, the points for which a large deviation of 
each ratio from 1 is seen are also generally the ones giving a relatively good y^ fit to the data. 

A similar trend is seen also for other values of However, deviations of Di and D 2 from 
unity by up to 40 — 50% are obtained for = 3° (fig. El^b)) and = 10° (fig. El^c)), but there 
are many more points with significantly large deviations of each of the ratios for the latter phase 
compared to the former one. For = 30° all the points appearing in fig. [2]^d) give T^i, D 2 and 
7^3 smaller than 1 and the overall deviation is generally smaller than for other non-zero phases but 
larger than for the rNMSSM limit. Thus, for this phase, the measured signal strengths can provide 
a clear indication whenever two Higgs bosons are present near 125 GeV instead of one. The reason 
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Figure 2: The ratios Di, D 2 and Ds, defined in eq. m, as functions of the difference between H 2 
and Hi masses, in the scenario when i^obs — Hi-\- H 2 . In (a) is set to 0°, in (b) to 3°, in (c) to 
10° and in (d) to 30°. The heat map in all the panels corresponds to the total • 
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why the deviations of the three ratios are much smaller overall in the case of (p,^ = 0 ° than for the 
CPV cases, for points showing the highest consistency with the data, will become clearer below. 

As noted earlier, a scenario with two mass-degenerate Higgs bosons in the cNMSSM entails both 
the i^obs — Hi -\- H 2 and i^obs = Hi/H 2 -\- Ai possibilities of the rNMSSM. Thus it is interesting 
to see which one of these two possibilities is favoured more by the data, for a given In the left 
panels of fig. [3] we thus show the squared normalised coupling Cyy against Cyy^ with the heat map 
corresponding to the total Similarly, in the right panels we have plotted Cyy vs. Cyy^ while 
the distribution of ttihs is shown by the heat map. For clarity of observation, we have included 
in this figure points with a total reaching up to 80, which is much higher than for the points 
shown in the earlier figures for this scenario. Also we have imposed an upper cut-off of 300 GeV on 
the mass of H^. We expect Cyi, to either vanish when a given Hi is a pure pseudoscalar (in the 
rNMSSM limit) or be relatively small when it is pseudoscalar-like (for > 0°). Note that these 
couplings satisfy the sum rule [26] 

N 

^VV — ^ 

i=l 

where N is the total number of neutral Higgs bosons that have a tree-level coupling to the gauge 
bosons, i.e., 5 in the cNMSSM and 3 in the rNMSSM limit H In the figure we see the above sum rule 
being satisfied almost completely by the three lightest neutral Higgs bosons under consideration 
here, implying that the remaining two doublet-like Higgs bosons are nearly decoupled. 

In the case of = 0° (i.e., in the rNMSSM limit) in the left panel of fig. El^a), we see two 
distinct kinds of points. There are some points lying along the diagonal, for which Cyy and Cyy 
alone are enough to satisfy the sum rule in eq. (|T^ . It is further evident from the right panel that 
Cyy for these points is exactly 0 . Hi and H 2 in these points should thus be scalars and H^ a 
pseudoscalar (i.e., Ai). But for the majority of the points, which lies along either of the axes, Cyy 
is nearly 1 , implying it is an almost pure doublet-like scalar, while Cyy is exactly 0 , implying it is 
a pseudoscalar, or vice versa. One can then observe in the right panel that for such points Cyy^ 
with Hs being the singlet-like scalar, is responsible for the sum rule being satisfied. Thus when 
the doublet-like scalar, whether Hi or i 72 , has Cyi^ slightly below 1, Cyy is slightly above 0. The 
mixing of the doublet scalar with H^ increases as its mass decreases, as is evident from the heat 
map in the right panel of the figure. As a result, the largest Cyy^ ^ 0 . 8 , is seen for the lowest tuhs 
obtained, which lies just above the allowed i^obs mass window. 

A closer inspection of the heat map in the left panel of fig. [3l^a) reveals that the lowest values 
of obtained for points lying along one of the axes, i.e., when the doublet-like scalar is nearly 

mass degenerate with the pseudoscalar. For points along the diagonal, the x^ is in fact always 
larger than 71. This is the reason for the relatively small deviations of D 2 and 7^3 from 1 seen 
in fig. El^a), where only points with x^ lower than 71 were shown. For such points, since one of the 
Hi and H 2 is a pure pseudoscalar as well as singlet-dominated, its contribution to the combined 
signal strength in the WW channel is null and that in the 77 and bb channels is minimal. Therefore, 
while the presence of Hi and H 2 of the rNMSSM near 125 GeV may possibly cause D 2 and 
to deviate more significantly from 1 , the consistency of this scenario with the LHC data is worse 
than that of the Hi -h Ai scenario. 

Fig. El^b) shows that, for = 3°, Hi and H 2 are almost always scalar-like while H^ is highly 
pseudocalar-like with a relatively much smaller Cyy generally. However, due to CP-mixing, Cyy 

^Note that since the /ism is a hypothetical SM Higgs boson with the same mass as a given Hi, at the tree level 
the ratio in fact corresponds to {gHiX/ghsM^)"^ ^^3 the equality in eq. ([T 3 l) is exact. However, since have 
actually been defined here in terms of the partial decay widths of Hi in the X channel, which include higher order 
effects also, the sum of may deviate slightly from unity. 


15 







0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Ot/ 


80 

78 

76 ^ 
o 

CO 

- 'I 

72 

70 



rn 300 
290 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

pH 


270 > 
a; 

260 

250 g 

240 

230 

220 


(d) 


Figure 3: Squared normalised coupling of Hi to the gauge bosons vs. that of H 2 (left) and of 
(right) in the scenario when i^obs = + ^2, with (j),^ set to (a) 0°, (b) 3°, (c) 10° and (d) 30°. 

The heat map in the left panels shows the distribution of in the right panels that of ttih^ . 
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can reach as high as 0.7 or so when the mass of is close to that of Hi and 7^2, though this 
happens for only a few points. A very crucial point to note here is that the total fhe left panel 
never falls below 68, which is due to the cut-off on the allowed upper value of This means 

that the points which give the overall best fit to the data have a much higher mass, which leads 
to a much smaller scalar-pseudoscalar mixing and hence negligible Cyy. 

For the = 10° case, illustrated in fig. El^c), while the maximum Cyy obtained is relatively 
small and hence Cyy and Cyy do not deviate from the diagonal by much in the left panel, there 
are many more points, compared to the (j),^ — 3° case above, for which Cyy is significant, according 
to the right panel. Finally, for = 30°, although Cyy never completely vanishes, it also stays 
smaller overall than it is for other phases. The reason for this is that the pseudoscalar-like 
never achieves a mass below 220 GeV or so, as can be noted from the heat map in the right panel 
of fig. [3](d). In the left panel one therefore sees that Cyy and Cyy always remain very close to 
the diagonal. Hence, for non-zero the data clearly favours two scalar-like Higgs bosons near 
125 GeV, instead of a pair of scalar-like and pseudocalar-like Higgs bosons. 

5 Conclusions 

In summary, we have tested the consistency of the real and complex NMSSM with the latest 
Higgs boson data from the LHC Run-I and the Tevatron. In particular, we have focused on 
scenarios wherein the resonant peak seen by the experiments can be explained in terms of two 
nearly mass-degenerate Higgs states around 125 GeV. Such scenarios have been verified in the 
rNMSSM previously and have not been ruled out yet. What we have shown here is that the 
possibility of such dynamics being available in the NMSSM is somewhat enhanced if some degree 
of (explicit) CP-violation is allowed in the Higgs sector. This can be done by assuming one or more 
of the Higgs sector parameters to be complex. By choosing this parameter to be one can evade 
the fermion EDM measurements, which tightly constrain the other possibly complex parameters in 
the Higgs and soft SUSY sectors of the NMSSM. 

In order to achieve the above we have performed detailed numerical scans of the parameter 
space of the cNMSSM to obtain various possible configurations with ^ 125 GeV Higgs boson(s) 
that also give SM-like signal strengths. In these scans we set the phase of to five different values, 
0°, 3°, 10°, 30° and 60°. Through a comprehensive analysis of the points obtained from these 
scans, we have then established that certain parameter configurations yielding two Higgs bosons 
near 125 GeV are slightly more favoured by the current data compared to scenarios with a single 
^ 125 GeV Higgs boson. This statement is even stronger when the two Higgs bosons are CP-mixed 
states. For the case of (j)^ = 10° we thus obtained: i) the point with the global minimum 
ii) more points with lying within 4 units of the global minimum compared to all other 
scenarios and phases tested; iii) more points with larger deviations of the ratios Di, D 2 and 
from unity. 

While analysing the aforementioned scenario with two Higgs bosons near 125 GeV, we have made 
sure that their masses are close enough that these two states cannot be distinguished experimentally 
as separate particles. In doing so we have exploited the fact that the experimental measurements are 
currently unable to reconstruct Breit-Wigner resonances, given that the experimental resolution in 
all channels investigated in the Higgs analyses is significantly larger than the intrinsic Higgs boson 
widths involved (so that LHC data actually reproduce Gaussian shapes). However, (tree-level) 
interference and (1-loop) mixing effects become crucial and need to be accounted for when the 
(pole) mass difference between two Higgs states is comparable or smaller that their individual 
intrinsic width. While we have ignored such effects here for points where they can be relevant. 
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which however make up a very tiny fraction of all the good points from our scans, they are the 
subject of a dedicated separate study [12]. 

Finally, in our analysis we have used up-to-date sophisticated computational tools in which 
state-of-the-art theoretical calculations and/or experimental measurements have been implemented, 
so that the solidity of our results is assured. 
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